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Abstract 

Changes in the medium of biological cell nucleus under ion beam action is consid¬ 
ered as a possible cause of cell functioning disruption in the living body. As the most 
long-lived molecular product appeared in the cell after the passage of high energy ions, 
the hydrogen peroxide molecule is picked out. The possibility of the formation of stable 
complexes of hydrogen peroxide molecules with active sites of DNA nonspecific recogni¬ 
tion (phosphate groups of the double helix backbone) is studied, and the formation of 
stable DNA-peroxide complexes is considered. Due to the negative charge on the oxygen 
atoms of DNA phosphate group in solution the counterions that under natural conditions 
neutralize the double helix have been also taken into consideration. The complexes con¬ 
sisting of oxygen atoms of DNA phosphate group, H 2 O 2 and H 2 O molecules, and Na"'' 
counterion have been considered. Energy of the complexes have been determined based 
on the electrostatic and van der Waals interactions within the approach of atom-atom 
potential functions. The stability of various configurations of molecular complexes has 
been estimated. It has been found that hydrogen peroxide molecules can form the stable 
complexes with phosphate groups of DNA and counterions which are no less stable than 
the complexes with water molecules. It is shown that the formation of stable complexes 
of H202-Na^-P0j can be detected experimentally by the observation of specific DNA 
vibrations in the low-frequency Raman spectra. The interaction of H 2 O 2 molecule with 
phosphate group of the double helix backbone can block the processes of DNA biological 
functioning and induce the deactivation of the genetic apparatus of the cell. Thus, the 
new channel of high-energy ions action on living cell has been proposed. 


1 Introduction 

Invention of new instruments for cancer therapy is known to be one of the most important 
challenges of modern science. In spite of the increase of funding for the development of new 
anticancer drugs the illness rate is still very high [1]. In this regard, the elaboration of non 
pharmacological methods of cancer therapy are stimulated. The last decades, use of high energy 
ion beams has established itself as an effective treatment of cancer disease. This method is based 
on well-known Bragg effect [2H5]. Heavy ions (usually protons or accelerated to the 
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energies of hundreds of MeV, are targeted into tumor. Within organism tissue the ions lose 
their initial energy mainly at the end of their trajectory of motion, and the energetic peak is 
formed (the Bragg peak). The initial beam energy is chosen in such a way that the position 
of the Bragg peak coincide with the location of tumor in human body (it can be 5-10 cm from 
body surface). The energy which is lost in living tissues destroys cancer cells and consequently 
the whole tumor. The method of ion beam therapy is much more effective than X-ray therapy, 
because of local action in the Bragg peak, especially in the case of places of difficult access, 
such as human brain M. In spite of great practical interest and large number of ion therapy 
facilities already built, there is no clear understanding of the molecular mechanism of action of 
high energy ions. 

Starting with the early study [6] the action of radiation on living organisms is known to be 
related with the DNA damage. The breaking of DNA double helix affects the mechanism of 
storage and transfer of genetic information which leads to cell death. Under heavy ion irradia¬ 
tion the breaks of hydrogen bounds in nucleotide pairs occur, which cause DNA melting, single 
and double strand breaks in double helix M. Single strand breaks in DNA and melting can 
be repaired due to the complementary structure of double helix and cell reparation mechanisms. 
The double strand breaks in DNA lead to damage of macromolecule, which can not be repaired 
[7]. Thus the formation of large number of double strand breaks of double helix is necessary 
for deactivation and destruction of DNA molecule. 

The double strand breaks in DNA macromolecule occur due to the large quantity of sec¬ 
ondary electrons born in the medium under the action of ion beam [8l-[Tn(IT3] . Some molecules 
decay into atoms and radicals that can break chemical bonds in biological macromolecules. 
At the same time, the calculations show that the number of the double strand breaks, made 
by secondary electrons and radicals, is insufficient for cell death |T0]. Therefore, additional 
channels of ion beam action on living tissues should exist. 

To describe DNA damage the formation of the shockwave in the medium after passage 
of high energy ions is suggested mma. Within the framework of shockwave mechanism the 
destruction of DNA double helix may occur even in the case of macromolecules situated far from 
ion hit point in the cell. The molecular dynamics simulations shows that shockwave cause large 
number of double strand breaks in DNA double helix [13]. However, in spite of high efficiency 
of proposed mechanism, there are processes which can prevent DNA damage in the living cell, 
particularly, the histone proteins restraining the macromolecule in frames of chromatin. 

The deactivation of DNA macromolecule may be induced by changes of medium that occur 
after the passage of high energy ions. The goal of the present work is to consider the changes 
of medium in the Bragg peak region and to pick out new molecular compounds that can 
influence the structure and dynamics of DNA double helix. The analysis of changes in water 
medium induced by ion beams is carried out in second section and the most probable molecular 
compounds formed in the Bragg peak are analyzed. The third section is dedicated to the 
description of calculation method of interaction of appeared molecules and atoms of DNA double 
helix. The results of energy calculations are presented in fourth section and the most stable 
complexes are selected. In the fifth section the low-frequency vibrations of formed complexes 
are studied, and the specihc modes of the complexes are found in the Raman spectra of DNA. 
The possible action of medium changes on DNA biological functioning is discussed in the sixth 
section. 
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Figure 1: Products of water radiolysis after 10 ® s of ion passage. Bar graph is built using the 
Monte Carlo simulation data [ 20 ] , 


2 Medium change under ion beam radiation 

The ions of high energy passing through intracellular medium induce the process of decay of 
water molecules (water radiolysis), different chemical and nuclear reactions [mi2T] . The pro¬ 
cesses of nuclear fragmentation are showed be not signihcant for environment [T5l[T6] , whereas 
the fragmentation of water molecules change essentially the composition and properties of water 
medium [3fl5|[Tn]. In the process of water radiolysis the secondary electrons (e“), radicals (OH* 
and HO*), molecules (H 2 O 2 and H 2 O), and other products (H*, OH“, H 3 O+) are produced 

[IM9]. 

Quantity and composition of radiolysis products depend on absorbed radiation dose. At 
small doses the number of secondary electrons and radicals is the largest, while the increase 
of the dose energy leads to the increase of molecular products yields. Due to the high energy 
transfer the number of molecular products is expected to be the largest in the Bragg peak 
region. The Monte Carlo simulation data for water radiolysis [201121] show that in the Bragg 
peak the secondary electrons and hydroxyl radicals yields prevail at early stage of radiolysis 
(time 10“^^ s). After the time 10“® s the reactions OH*-|-OH* ^H 202 , HO*-|-HO* ^H 202 -l -02 
occur, and the yields of molecular products become larger than other yields IZD]. Taking into 
consideration that the microsecond time scale is typical for biological processes, the equilibrium 
composition of the intracellular medium after the passage of high energy ions must be enriched 
by the molecular products of water radiolysis. 

To analyze the composition of the radiolysis products at the Bragg peak the Monte Carlo 
simulation data [20] are used. The distribution of radiolysis products at the time 1 //s after ion 
passage are showed in the Figure 1. It is seen that the number of hydrogen peroxide molecules 
in intracellular medium is larger than other products of water radiolysis. The number of H 2 
molecules is also considerable, but these molecules don’t have high reactivity and as a result do 
not threaten biological macromolecules. Radicals which also have substantial yields after some 
time transform into neutral water molecules or hydrogen peroxide. Thus, the composition of 
intracellular medium changes after ion beam irradiation and signihcant amounts of hydrogen 
peroxide molecules appear. 
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Table 1: Parameters of potentials. 


Atomic charges in molecules 

Potential (4) [29] 

Potential (3) 

[2411^ 


Molecule 

H 0 

Parameter Na’*' 

Parameter 

0-0 

0-H 

H2O 

H2O2 

PO4 

+0.33 (e) -0.66 (e) 
+0.41 (e) -0.41 (e) 
-0.50 (e) 

q (e) +1.00 

ro (A) 2.35 

b (A) 0.3 

Aij (A*^kcal/mol) 
Bij (A^^kcal/mol) 

200 

12900 

86 

31300 


The concentration of hydrogen peroxide in the living cell is kept constant by special fer¬ 
ments [22], but it may be increased greatly by the action of ion beams. Increasing the amount 
of hydrogen peroxide molecules raises the probability of their interaction with DNA. The inter¬ 
action of large number of H 2 O 2 molecules and DNA can deactivate the cell genetical apparatus 
and induce cell death. Therapeutic effect of the hydrogen peroxide has been found under the 
treatment of cancer disease (see for example |23|)- Ion therapy allows to increase the concen¬ 
tration of hydrogen peroxide in tumor cells with high precision and with minimal damage of 
healthy tissues. 

The hydrogen peroxide molecule can interact with DNA molecule via charged atomic groups 
of the double helix. It is well known that the phosphate groups of double helix backbone bear 
the highest charge [7]. Under the natural conditions it is neutralized by metal ions (for example, 
Na"*" or K+). Thus, the interaction of H 2 O 2 molecule with the phosphate groups of the double 
helix and tethered metal counterions is expected to be the most probable. 


3 Model and method of energy calculation 

The capability of hydrogen peroxide to interact with DNA is studied by calculating the energy 
of different complexes. The considered complexes consist of DNA phosphate group, hydrogen 
peroxide molecule, water molecule and Na^ counterion. For the determination of relative 
stability of different complexes two oxygen atoms of the phosphate group are sufficient to 
consider. Other atoms of POj are essential for accurate calculations of energy. Schematic 
constructions and structural parameters of molecules in complexes are showed at Figure 2. 

The method of atom-atom potential functions is used for calculation of energies of complexes 
[211125]. Within the framework of this method the energy is presented as the sum of pair 
interactions between atoms which belong to different structure elements. The energy of complex 




Figure 2: Schematic structure of hydrogen peroxide molecule, water molecule, and phosphate. 
Characteristic values of bond lengths, valence angle and dihedral angles are showed. 
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can be presented as the snm of two components: 

U = [Ueijrjj) + Uvdwjrij)] , ( 1 ) 

ij 

where Uei{rij) and U^dwirij) are the potential energy of electrostatic and van der Waals inter¬ 
actions between i and j atoms at the distance Vij. 

Electrostatic interaction is dehned as follows: 


Uel{Tij) 


QiQj 

4:7reeorij ’ 


( 2 ) 


where g* and qj are the charges of i and j atoms, Eq is the dielectric constant of vacnnm, e is the 
dielectric constant of medinm. Charges of atoms of hydrogen peroxide and water molecnles are 
determined nsing the valnes of dipole moments 2,1 D (H 2 O 2 ) and 1,86 D (H 2 O) |26]. Charges 
of atoms in phosphate gronp are taken only for oxygen atoms in order to eqnate total charge 
with -e (Table 1). The valne of dielectric constant e = 1 is nsed in calcnlations. 

In the case of van der Waals interaction the Lennard-Jones potential is nsed: 


U^dw i.^ij ) 


^-6 




-12 


(3) 


where Aij and Bij are the constants (Table 1). The hrst and the second terms in (3) describe 
the van der Waals repnlsion and short range repnlsion between atoms, respectively. To take 
into consideration the hydrogen bonds (0-H- • -O) the interaction between respective oxygen 
atoms are calcnlated using the modihed potential [25], where exponent in the hrst term of 
equation (3) is changed to -10. Such exponent describes the short-range character of hydrogen 
bonds. 

For the description of interaction between ion and charged oxygen atoms the Born-Mayer 
potential is used: 


UsMir ij) 


mj 

iTTSEorij 



(4) 


where b is the repulsion constant, tq is the equilibrium length (Table 1). This potential was 
developed to describe the ion crystal energy, and it takes into account electrostatic attraction 
and short-range repulsion The potential (4) has been adapted in [281129] for the description 
of counterion interactions with the phosphate groups of DNA double helix backbone. 

Using the formulae (1) - (4) with parameters in the Table 1 the energies of different com¬ 
plexes may be calculated. 


4 Complexes of hydrogen peroxide with DNA 

Within the framework of proposed method the energies of two-component (H 2 O 2 -PO 4 , H 2 O- 
PO 4 , Na’''-P 04 , H 202 -Na’'', H 20 -Na''') and three-component (H 202 -Na’''-P 04 , H 20 -Na'''- 
PO 4 ) complexes are calculated. The most stable conhgurations for each complex are found. 

Two-component complexes. The most stable state for the complex Na^’^-TOj is obtained in 
case when sodium ion is localized at the distance 2 A from the center of 0-0 segment of the 
phosphate group. The calculated anergy is the lowest among two-component complexes (Fig. 
3). This complex corresponds to neutralization of DNA phosphate group by counterion under 
the natural conditions. 
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-12 kcal/mol 

Figure 3: The structures of two-component complexes: 
H2O2-PO4, H2O-PO4 . 



11 kcal/mol 



Na+-P 04 , H 202 -Na+, 


H20-Na+, 


The most stable complex of sodium ion with hydrogen peroxide is obtained in the case of 
ion localization at distance 2.66 A to the center of 0-0 bond with equal distances to H atoms 
on the opposite side of the molecule (Fig. 3). In the case of Na’''-H 20 the most stable complex 
is obtained for the Na’''-0 distance 2.74 A. Under the natural conditions water molecules in 
complex with Na^ are the part of ion hydration shell. The hydrogen peroxide molecule may 
substitute the water in the ion hydration shell. 

The equilibrium configurations of H 2 O 2 -PO 7 complex is obtained in the case of symmetrical 
localization of H atoms of hydrogen peroxide molecule with respect to the 0-0 segment in the 
phosphate group (Fig. 3). The complex has minimal energy for the distance 2.85 A between 
centers of 0-0 segments in the molecule and in the phosphate group. In the case of H 2 O-PO 7 
complex the equilibrium conhguration is obtained for the distance 2.9 A of water oxygen to 
0-0 segment in the phosphate group (Fig. 3). 

The hydrogen bonds, which can be formed between H (water or hydrogen peroxide) and O 
(phosphate group), do not influence signihcantly the energy values. In the case of the strongest 
hydrogen bond (O-H-0 atoms are localized on the line) the calculated energies are lower 
than energies of Il 202 -Na+-P 04 and Il 20 -Na+-P 04 complexes. In configurations with two 
hydrogen bonds between atoms of II 2 O 2 ( 112 ©) molecule the O-H-O segment is bent more than 
30°, that makes such hydrogen bonds very weak. 

The equilibrium configurations of II 2 O 2 -PO 7 and II2O-PO4 complexes depend on mutual 
disposition of molecules. Rotating H 2 O 2 molecule with respect to the axis passing through the 
center of 0-0 segment, the most stable complex is found for the angle about 60° (energy - 12.8 
kcal/mol) (Fig. 4a). In the case of molecule rotates around 0-0 bond, the energy minimum is 
obtained for the angle 0°, while at 180° the repulsion occurs (Fig. 4b). Some small minimums 
at about the angles 70° and 290° appear due to the van der Waals repulsion and Coulomb 
attraction (dash lines at figure). In the case of rotation of water molecule the obtained results 
are qualitatively similar. 

Three-component complexes. The calculations are carried out for Il 202 -Na^-P 07 and 
H 20 -Na’''-P 07 complexes. The energy of the complexes are calculated for different distances 
between hydrogen peroxide (water) molecule and ion-phosphate complex. The ion-phosphate 
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Figure 4 : Energy of H2O2-POJ complex as a function of rotation angle of H2O2 molecule, 
a) The dependence of energy on rotation angle with respect to the axis passing through the 
centers of 0-0 segments in H2O2 and PO4 (is showed in inset), b) The dependence of energy 
on rotation angle with respect to the 0-0 bond of H2O2 molecule (is showed in inset). The 
van der Waals and Colomb contributions to the energy are showed by dashed lines. 


distance is taken equal to 2 A, correspondingly to the result for two-component complex Na"*"- 
PO4 (Fig. 3 ). Such ion localization is the most energetically favorable, therefore other variants 
of ion localization are not considered in the present work. 

The conhgurations where hydrogen atoms of H2O2 (H2O) molecule are turned toward and 
outward phosphate group are considered (Fig. 5 a). In the hrst case there is no energy minimum 
due to the electrostatic repulsion between hydrogen atoms and sodium ion. Increasing the 
distance to H2O2 (II2O) the energy of the complex monotonously decreases to the value of 
Na’''-P04 complex (Fig. 5 b). In the case when H atoms of H2O2 (II2O) molecule are turned 
toward the phosphate group the energy minimum 130.3 kcal/mol ( 129.2 kcal/mol) appear at the 
distance 2.75 A ( 2.82 A). The minimum is conditioned by the electrostatic attraction between 
ion and oxygen atoms of molecule. Rotation of hydrogen peroxide or water molecules around 
axis which passes through center of 0-0 segment plays insignihcant role because the total 
energy depends on interaction between ion and oxygen atoms mainly. Thus, the stability of 
H202-Na+-P04 and H20-Na’''-P04 complexes are almost equal. 

5 Vibrations of hydrogen peroxide in DNA low-frequency 
Raman spectra 

The formation of stable complexes of II2O2 molecule with DNA phosphate groups should be 
accompanied by the appearance of specihc modes of vibrations of hydrogen peroxide molecule in 
the vibrational spectra. The frequencies of vibrations of Il202-Na+-P07 complex are expected 
to be in the low-frequency spectra range (<200 cm“^), analogically to the modes of counterion 
vibrations with respect to DNA phosphate groups (ion-phosphate modes) [ 28 II 3 T] . The low- 
frequency spectra of DNA is known to be characterized by the conformational vibrations of the 
double helix structural elements as whole (nucleosides, nucleotides, phosphate groups) [ 32 II 35 ] . 
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Figure 5 : Energy of three-component complexes, a) Position of H2O2 and H2O molecules in 
the considered complexes, b) Energy of the complex as a function of distance between H2O2 
or H2O molecules to the ion-phosphate complex. 
































The vibrations of hydrogen peroxide molecnle with respect to the DNA phosphate gronps shonld 
distnrb the internal dynamics of the donble helix changing the low-freqnency spectra shape. 

To determine the inflnence of H 2 O 2 molecnles on DNA conformational vibrations the phe¬ 
nomenological model describing the Raman low-freqnency spectra of DNA with connteri- 
ons [ 2 HH 3 O] is used. Within the framework of this model the DNA double helix is repre¬ 
sented as the double chain of nucleotides with counterions tethered to the phosphate groups 
of the macromolecule backbone. To take into consideration the influence of hydrogen peroxide 
molecule this model is modified by increasing the mass of tethered counterions for the mass 
of one H 2 O 2 molecule. The low-frequency vibration spectra of DNA with counterions is also 
calculated for the case of H 2 O molecule attached to the ion. The equations, necessary for the 
calculations of the mode frequencies and Raman intensities, are the same as in the original 
works [29l[3n] . The calculation parameters are taken for the case of R-form of the double helix. 

As the result the low-frequency Raman spectra of DNA with attached H202-Na’'' complex 
is calculated as well as the spectra of DNA with H20-Na’'' complex (Fig. 6). The internal 
vibrations of the double helix (lower than 120 cm“^) are connected with the vibrations of H- 
bond stretching in nucleic bases and the vibrations in nucleosides due to the flexibility of sugar 
ring. These modes do not sense essentially the presence of hydrogen peroxide molecule. At the 
higher frequencies the modes, characterizing the vibrations of the complexes with respect to 
the phosphate groups, are determined about 155 cm“^ and 140 cm“^, respectively. The mode 
intensity is higher in the case of complex with hydrogen peroxide molecule. Thus, our estimates 
show that hydrogen peroxide molecule, attached to DNA phosphate groups, may induce the 
visible changes in the low-frequency Raman spectra. 





Figure 6: The low-frequency Raman spectra of DNA with H202-Na’'' complex (solid line) and 
with H20-Na’'' complex (dashed line). In the insets the structure of respective complexes are 
showed. 
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6 Discussion 


The calculations show that the complexes including hydrogen peroxide and water molecules 
are the most stable under the presence of sodium counterion. Without Na"*" the stability of 
the complexes decreases for about one order. That is due to the dominant role of electrostatic 
interactions in the complex formation. The stability of H202-Na+-P0j is about the same 
as H20-Na+-P0j. Therefore, hydrogen peroxide molecule may be expected to incorporate 
into the hydration shell of sodium counterion already tethered to the oxygen atoms of DNA 
phosphate group. 

The formation of stable complex of H202-Na’'' with the phosphate group may be detected 
by the low-frequency Raman spectra of DNA. According to our estimations the vibrations of 
counterion with hydrogen peroxide molecule with respect to the phosphate group of the double 
helix backbone should be about 10 cm“^ lower than the frequencies of vibrations of the complex 
with water molecule. Taking into consideration the determined frequency shift the lifetimes of 
H202-Na+-P0j and H20-Na+-P0j complexes may be compared. According to the Arrhenius 
equation the ratio between characteristic lifetimes may be written in the following form: 


Ti a;2 f El — E2\ UJ2 

— = — exp —^ 

T2 Ui \ kl J Ui 


(5) 


where ri, Ui and r 2 , UJ 2 are the lifetimes and frequencies of the complex with water and hydrogen 
peroxide molecules, respectively; Ei and E 2 are the energies of the complexes with hydrogen 
peroxide and water molecules, respectively. The energies of the considered complexes are rather 
close (Fig. 5), therefore the value of the energy difference can be equal to zero. As the result the 
expected lifetime of the complex with H 2 O 2 is higher than in the complex with H 2 O molecule. 

Due to the larger lifetime of H 202 -Na’''-P 04 complex the hydrogen peroxide molecules may 
be accumulated near DNA double helix. At signihcant amounts the H 2 O 2 molecules can influ¬ 
ence DNA biological functioning by two mechanisms. Firstly H 2 O 2 can induce double strand 
breaks in double helix backbone by decaying into radicals, and secondly H 2 O 2 molecules may 
break nucleic-protein recognition, stopping the processes of translation of genetic information. 
Thus, the hydrogen peroxide molecule, interacting with DNA double helix, can deactivate cell 
genetic apparatus and induce cell death. 


7 Conclusions 

The hydrogen peroxide molecules appear in the cell medium after the passage of high-energy 
ions. To study the possible effect of H 2 O 2 molecule in the functioning of the living cell the 
complexes of hydrogen peroxide molecules with the phosphate groups of DNA backbone has 
been considered. Using the method of atom-atom potential functions the stability of different 
complexes, consisting of DNA phosphate group, H 2 O 2 (or H 2 O) molecule, and Na"*" metal ion, 
have been studied. The results show that the complexes with H 2 O 2 molecule and phosphate 
group of DNA are no less stable than respective complexes with H 2 O molecule. The sodium 
counterion neutralizing the charge of phosphate group plays the key role in stabilization of these 
complexes. For the case of most stable complexes of hydrogen peroxide with the phosphate 
groups the low-frequency Raman spectra of DNA have been calculated. The vibrations of 
the complex as whole with respect to DNA phosphate group has been found near 150 cm“^. 
The observation of this mode may be considered as a test revealing the complex formation. 
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Taking into consideration the determined frequencies of vibrations the lifetimes of H202-Na’''- 
POj and H 20 -Na’''-P 0 j complexes have been compared, and the lifetime of the complex with 
hydrogen peroxide has been found higher than in the case of water molecule. Taking this 
into consideration the hydrogen peroxide molecules are expected to be accumulated near DNA 
macro molecule. The action of H 2 O 2 on DNA structure and dynamics may be sufficient for 
the braking the processes of DNA biological functioning. The influence of hydrogen peroxide 
molecules on DNA can be regarded as an individual channel in heavy ion therapy. 
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